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DETAILED COMPUTATIONAL PROCEDURE FOR DESIGN OF CASCADE BLADES WITH PRE- 
SCRIBED VELOCITY DISTRIBUTIONS IN COMPRESSIBLE POTENTIAL FLOWS 1 

By George R. Costello, Robert L. Cummings, and John T. Sinnette, Jr. 


SUMMARY 


SYMBOLS 


A detailed step-by-step computational outline is presented for 
the design of two-dimensional cascade blades having a prescribed 
velocity distribution on the blade in a potential flow of the usual 
compressible fluid. This outline is based on the assumption 
that the magnitude of the velocity in the flow of the usual com- 
pressible nonviscous fluid is proportional to the magnitude of 
the velocity in the flow of a compressible nonviscous fluid with 
linear pressur e-volume relation. The computational procedure 
includes several ways of adjusting the prescribed velocity to sat- 
isfy restrictions imposed by the method. Tables of coefficients 
are given for evaluating the necessary integrals, including the 
determination of the harmonic conjugate function. Numerical 
examples are included. 


INTRODUCTION 

A method for computing blade profiles with prescribed 
velocity distributions based on the assumption that the pres- 
sure-volume relation is linear is presented in reference 1. 
The method uses the prescribed velocity distribution and 
compatible free-stream conditions to determine a mapping 
function that transforms an incompressible flow about the 
unit circle into an exact compressible flow, with linear pres- 
sure-volume relation, about a cascade of blades having the 
desired velocity distribution. 

In order to apply this method to the design of a cascade 
with a given velocity distribution and free-stream conditions 
in a flow of an “actual” fluid (the usual compressible non- 
viscous fluid), flic relation between the actual fluid and the 
fluid with the linear pressure-volume relation must be approx- 
imated so that the required velocity distribution and free- 
stream conditions for the second fluid may be determined. 
In this investigation, which was made at the NACA Lewis 
laboratory, the magnitudes of the velocities in the two fluids 
are assumed to be proportional and the constant of propor- 
tionality is determined b_y tlie continuity equation using the 
same upstream and downstream flow angles for the two fluids. 

By use of this relation, a computational procedure was 
developed to obtain the blade profile with a minimum of 
effort. The procedure, presented herein, includes the adjust- 
ment of the prescribed velocity distribution to satisfy the 
restrictions oil the mapping function and the numerical 
computation of the harmonic conjugate function. 


The following symbols are used in this report : 

A, B, C 3 , . . . constants 

a ,d 

a, b, c constants 

C(e ie ) function of d defined by equation (18) 

d spacing of cascade blades 

H(e ie ) function of 8 


H(e ie ) 

Im 

im 

im 

k 

M 

n 

Q(s) 

2 


~t 

Re 

A(ReH) 

r 

s 

s 

u 


v{.e) 

z=x+iy 

a 

V 

7 

S 

t 

e 


final values of H{e id ) 
imaginary part 

distorted velocity defined by equation (22) 
modified Kid) 

constant determined by range of potential 
Mach number 

constant determined by trailing-edge angle of 
blade 

dimensionless velocity on blade in direction of 
increasing s 

magnitude of dimensionless velocity (ratio of 
velocity to stagnation velocity of sound) for 
flow with linear pressure-volume relation 
final value of q after modification of ReIT(e ie ) 
real part 

change in ReH(e ie ) 

ratio of velocities (determined by equation (4)) 
arc length along blade measured counterclock- 
wise from tail (figs. 1 and 2) 
arc length along final blade shape 
magnitude of dimensionless velocity (ratio of 
velocity to stagnation velocity of sound) 
for actual fluid 

velocity on unit circle (incompressible flow) 
complex variable in cascade plane 
angle of velocity in compressible flow (meas- 
ured from positive z-axis) 
circulation (positive counterclockwise) 
ratio of specific heats 

difference between Re H(e i0 ) and parabola at 
half-interval point 

including trailing-edge angle of blade 
circle angle (incompressible flow plane) pos- 
itive counterclockwise 


• Supersedes NACA TN 2281, “Detailed Computational Procedure for Design of Cascade Blades With Prescribed Velocity Distributions in Compressible Potential Flows’’ by George R. 
Costello, Robert L. Cummings, and John T. Sinnette, Jr., 1951. 
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X auxiliary function defined by equation (15) 

p density 

r variable of integration 

ip velocity potential 


Subscripts: 

c 

1 
n 
t 
1 

2 


compressible flow 
incompressible flow 
leading edge (nose) 
trailing edge (tail) 
upstream 
downstream 


RELATION BETWEEN FLOWS 


In the present problem of designing a cascade, the pre- 
scribed conditions are the upstream velocity Uie ia i, the 
downstream velocity u 2 e ia a, and the magnitude of the di- 
mensionless velocity on the blade as a function of the arc 
length u=u(s) (see figs. 1 and 2). The upstream and down- 
stream velocities are related by the continuity equation 


A/ose 



Figure 2. — Prescribed velocity distribution on cascade blade. 


continuity relation for the flow with the linear pressure- 
volume relation: 

— - ■ ■ ■ <]i COS «1= COS Ot 2 (3) 

a 1 + <7r’ a 1 + '7i i 


( 



-u, COS «! = 


1 


lizl 

2 



u-> cos a 2 


(I) 


Consequently, only three of the quantities u u u 2 , a,, and 
may be assigned, and the fourth is determined by equa- 
tion (1). 

In utilizing the method of reference 1 to design the cas- 
cade, the prescribed conditions are employed to determine 
the upstream and downstream velocities and the velocity 
distribution on the blade for the fluid with the linear 
pressure-volume relation by assuming that the magnitudes 
of the velocities for the two flows are proportional; that is, 

q=ru (2) 

The constant of proportionality is determined from the 
additional assumption that the flows have the same free- 
st ream directions by substitution of equation (2) in the 


, cos 2 <*■> — Ui 2 cos 2 ai ,, . 

? =-y 77 a ;; — c (4) 

ay ad (cos- a, — cos- a 2 ) 

Hence, the upstream velocity qie' a i, the downstream vel- 
ocity <7 2 e‘“ 2 , and the velocity distribution on the blade 
g=qr(s) are now determined for the flow with the linear 
pressure-volume relation. 

Other approx-imations to the relation between the two 
flows could be used but, with the present cascade inter- 
pretation, the upstream and downstream flow angles and the 
spacing are the same for the two flows. The two density- 
velocity relations are also in good agreement, as shown in 
figure 3 where the density-velocity relations are plotted for 
Ui = 0.700 and u 2 = 0.500. When (he upstream and down- 
stream velocities are equal, the density-velocity relation, as 
given by this interpretation, is the same as the density- 
velocity relation in the Karman-Tsien approximation. 




Figure 3.— Comparison of donsity-velocity relations for mi= 0. 700 and «2=0. ">00. 
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NUMERICAL PROCEDURE 

After the free-stream conditions and velocity distribution 
on the blade are determined for the flow with the linear 
pressure-volume relation, the method of reference 1 is used 
to obtain the blade shape. In this method, the given 
conditions are used to select an incompressible flow about 
the unit circle and to determine the mapping function that 
transforms this flow into a compressible flow with the circle 
transforming into the desired blade. A detailed procedure 
for determining the circle flow and the mapping are given in 
the following paragraphs. 

FLOW IN CIRCLE PLANE 

The velocity potential < p c (s) and the circulation 1\ on the 
blade are obtained by integrating the prescribed velocity 
distribution q(s). For this purpose, it is convenient to 
define Q(s) by 

Q(/j)=z — q(s) 0 < <s < .s „ 

Q(s) = < ?(*) «„<* < 2 tt 

then 

= f QA) < I x (5) 

./<> 

r c = P' q ( n) ,/« (6) 

./o 



Figure 4.— Velocity and potential distribution for fluid with linear pressure- volume relation. 


where the total arc length has been taken equal to 2w. 
These integrals are evaluated by Simpson’s rule using suf- 
ficient points so as to obtain accurate values of <p c (s). The 
values of <p c (s) should be plotted on the graph with Q(s) 
(fig. 4) in order to obtain q easily as a function of the circle 
angle later in the computation. The spacing (1 of the 
cascade is given by 

d=- 7 5s 7 (7) 

2i sm a, — q 2 sin a 2 

From equation (35) of reference 1, the incompressible 
potential function for points on the unit circle may be 
written 

<Pi(d )=- 2 Re /ltanh -1 ^^ + (/m A+ImB) tan " 1 + 

(I m A— I in B) fan - 1 + D (8) 

where D is chosen to make the potential zero at the tail 

, . ir sin 6 

stagnation point d,, the angle convention is— ^<tan ‘ gj — ^ 

<-» and tan -1 r — 1 is taken in the same quadrant and the 
2 tan 1 1 k 

same direction as 6. The values of Re A, Im A, Re B, and 
hn B are determined from the free-stream conditions and 
spacing by 


He A ' ? ' rZ / |c ° S “’ 1 

2ir-\ 1 +2r 

(9) 

Im ^ 1 = — \ 1 + t?i 2 Re A tan a. 

(10) 

Cb 

II 

1 

Os 

(11) 

I m B = a 1 + ( h 2 He A tan a 2 

(12) 

The value of the constant k is determined by the condition 
that the potential range on the circle is equal to the potential 
range on the blade: 

<p c (27r) — <ft,( 6 '») = <Pi(0t + 2ir) — iPiid,,) 

(13) 

The proper value of k is computed by assuming a 
and computing the stagnation points d, and d„, 
the roots of the equation 

value for k 
which are 

. , , . — ( Im A+Im B) sinh k 
sm (0 + X) 2 Re A sec X 

(14) 

where X is given by 


• 

, (Im A—Im B) tanh k 
X = tan - 2 R ~- A 

(15) 

-1< X <I 
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and then determining the value of <p i (d t -\-2ir) — <p i (d n ). This 
process is repeated for different assumed values of k to ob- 
tain aplot of <p i (0<+2'!r) — <Pi(d n ) as a function of k. The value of 
k satisfying equation (13) is then obtained by interpolation. 
A computational outline of this process for use with a desk 
calculator is given in table I, in which the usual column op- 
erations are written as lines. 

When k has been determined, the flow about the circle 
is known with the value of the potential on the circle given 
by equation (8) and the velocity on the circle given by 


where 

(7 3 =ln ~^ 1 ■\lRe 2 A J r Im 2 Ae k [e~ k ( 2 cosh k — 

' <Zi 

2 cos 0()]~” /2 j 

C t = In j 1 + a- 1 ±g jf y >Re*B+Im 2 Be k [e-' : (2 cosh k + 

1 <Z-2 

2 cos 0 ( )] _,!/2 | 


4 Re A cosh k sec XT . . . , . . . ,1 ,, 

SS2J-C5T2-#- L“‘ + < 18 > 

A computational outline for v{8) and <p t {8) based on 80 points 
is given in table II, lines 1 to 19. 


(? 6 =tan 1 


C 6 =tan _1 



+ ai — n tan -1 1 
4 -a -i—n tan -1 


— sin 6, \ 
e k — cos 0, / 

' _ sinfl, \ 

y* + cos 8,J 


+ ir 


DETERMINATION OF MAPPING FUNCTION 

The mapping function (equation (62), reference 1) may be 
written in the form 


2 = 



v{8) 2 { 2 cosh 2k — 2 cos 2 g) i(e+|+«-tan-> 
e k -\j2 cosh 2k — 2 cos 28 


where 


si n 26 \ 
cos 2 9-e- k )(j q 

(17) 


g ]= _J 2 2 COS ( d t e)] »/2 e [Re C(e -»)+(2 cosh 21—2 cos 29) Re F/(«**)] 

g 2 —n — -'j+Im C(e itr ) + (2 cosh2& — 2 cos2 8)ImH(e ie ) 

and 

71=1-- 


where e is the trailing-edgc angle of the blade in radians. 
The mapping function is completely determined when 
Re C(e ie ), Im C(e i9 ), Re H(e ie ), and Im H(e ie ) are known. 

Determination of Re C(e ie ) and Im C(e ie ). — The function 
C{e tB ) is defined by equation (44) (reference 1) and may be 
written in the form 


<7(e i9 )=.I( l +«*-«) In 
i (1 — e k ~ ie ) In — 


(1 + 

-Be ia A k (i + Vl + g 2 2 ) 
(1 -\-e ie t ~ k ) n q 2 


(18) 


Hence Re C(e ie ) and Im C(e is ) are obtained by taking the real 
and imaginary parts of equation (18): 


/M?(e i9 )=i[(0 3 + ^) +«*(£,- Q cos e + e k (a-C 6 ) sin 0 ] - 

(19) 

ImO(c i9 )=i£(C 5 + C 6 ) + e*(C 5 — C 6 ) cos 8—e k (C 3 —C 4 ) sin 0 J 

( 20 ) 


The quadrants for the arc tangent terms in C 5 and C s are se- 
lected by considering the numerator and the denominator of 
the argument as signed quantities and choosing a quadrant 
in which the sine has the same sign as the numerator and 
the cosine has the same sign as the denominator. For con- 
venience, these angles are taken as positive and the resulting 
values of C b and C 6 are changed by multiples of 2i r until |C 5 | 
and | (7 6 [ are each less than 2ir. 

Determination of Re H (e i9 ). — In order to obtain Re II (e te ), 
the prescribed velocity along the blade must be determined 
as a function of the circle angle q=q(8). This relation is 
obtained from the equality of the potentials at corresponding 
points in the two planes. Because Q(s) and <p c (s ) are plotted 
on the same graph (fig. 4), the magnitude of the prescribed 
velocity at a given circle angle 8 may be found directly by 
reading the absolute value of Q (<Z= IQI) at the abscissa 
where <p c (s ) equals the calculated value of <pt(6) (table II, 
line 19). The prescribed velocity is thus obtained as a 
function of the circle angle q=q{8). 

From equation (48) (reference 1) 


In 

ReH(e tl! ) = - 

0 (0) K 2 cosh 2lc— 2 cos 20 

-ReC(e te )-\-k 

_ K (8) [2 — 2 cos (8 t — 0)] n/ 2 _ 

2 cosh 2k — 2 cos 20 

(21) 


where 

m- 2q(e) . 
1 +1 1 +<z(0) _ 

(22) 


Lines 20 to 36 of table II show the detailed computation for 
obtaining Re H(e ie ), including the evaluation of Re C(e ie ). 

Adjustment of Re H(e ie ).' — Restrictions on the mapping 
function require that Re II(e is ) satisfy the conditions 


J* ReH(e iS )d8= 0 

(23) 

j" Re H(e ie ) sin 0 dd=0 

(24) 

|* ReII{e ts ) cos 0 dd=0 

(25) 
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These integrals are evaluated numerically (in the manner 
developed by Dr. Glenn H. Peebles, while at Douglas Air- 
craft Co., Inc., in his work on isolated airfoils in incom- 
pressible flows) by using a parabolic variation of Re H(e ie ) 
between three consecutive points and integrating analyti- 
cally, : The , point spacing used in the examples of this 
report is an equal spacing of 80 points (designated whole 
points and numbered 0, 1, 2, ... 79) around the circle 
beginning with point 0 at 0= — 180° and continuing around 
the circle at 4.5° intervals. Where Re H(e‘°) has large 
fluctuations (usually in the neighborhood of the stagnation 
points), the midpoints of these intervals (called half-points) 
are also used. The coefficients for evaluating the integral in 
equation (23) are the Simpson’s rule coefficients. The 
coefficients for the integrals in equations (24) and (25), 
which are obtained by integrating analytically the product 
of the parabolas and sin 0 or cos 6, are given in tables III to 
VI. When half-interval spacing is used, it should begin and 
end at even-numbered points of the original 80 points and 
the integration coefficient to be used at the beginning (or 
ending) even-numbered point is one-half the sum of the 
whole-point coefficient and the half-point coefficient at that 
point. 

If the values of Re Tl(e le ) from equation (21) do not satisfy 
equations (23), (24), and (25), the velocity distribution 
chosen is incompatible with the chosen free-stream conditions 
and must be modified to make these integrals zero. The 
most desirable method of modification will depend on the 
magnitude of the integrals in equations (23), (24), and (25) 
and on the features of the original prescribed conditions that 
are to be preserved. Usually the original frcc-strcam con- 
ditions are retained and the velocity distribution is modified. 
When the integrals are large, however, changing the circle 
flow may be desirable in order to minimize the changes in 
the essential characteristics of the original velocity distribu- 
tion (such as limits on maximum velocity, diffusion rate, 
and so forth). 

Modification of the circle flow provides essentially one 
additional degree of freedom, as can readily be seen from 
equations (9), (10), (11), and (12). The strength of the 
singularities A and B are determined by the free-stream con- 
ditions and the spacing. The range of potential is deter- 
mined by the location of the singularities, that is, by the 
constant k. Only the ratio of the potential range to the 
spacing is important, however, because changing both to- 
gether merely changes the dimensions of the cascade by a 
scale factor. This additional degree of freedom may be con- 
veniently represented by the ratio of arc lengths of the upper 
and lower surfaces. Consequently, this ratio has a large in- 
fluence on the size of the integrals in equations (23), (24), 
and (25). The selection of the ratio may be based on the 
ratio for a blade having a similar velocity distribution; or, 
in some cases, it may be advantageous to try several ratios 
of arc lengths and roughly approximate the integrals using 
only 20 points and from these results select the proper ratio 
to minimize these integrals. The final adjustment to reduce 


these integrals to zero may then be obtained with a fixed 
circle flow by adjustment of Re H(e ie ). 

When adjustments are being made in Re H(e ie ), the change 
in profile arc length and profile velocity q(6) produced by 
these adjustments must be considered. The local arc length 
ds corresponding to dO will be decreased in the same ratio 
that the velocity is increased in order to maintain the same 
potential. The change in arc length is automatically ac- 
counted for in the final integration for the blade coordinates. 
The change in the velocity produced by changing Re H(e ie ) 
by an amount A (Re H) is indicated by the change in the 
“distorted” velocity K(d), which is related to q(8) by equa- 
tion (22) (fig, 5). The changed value, denoted by K(8), is 
given by 

(0) = _KT (0) g-(2coah 2fc— 2 cos 20) A (Re H ) (26) 

Making a change in Re H(e ie ) therefore has the effect of mul- 
tiplying the distorted velocity by a factor. This factor varies 
widely with the circle angle and produces the smallest change 
in K{6) at 0= ±t or 0=0 and the largest change at 0= ±ir/2. 
Hence, for a given percentage change in K{8), the greatest 
values of A(7?e II) occur near 0= or 0=0. The effect of 
k on these values of A {Re II) for a given ratio of K{d)/K(d) 
is shown in figure 6. 

Consequently, for small values of the integrals in equations 
(23), (24), and (25), complete adjustment can frequently be 
made by merely changing the slope of the distorted velocity 
K(d ) at the stagnation points. The changes in the integrals 
produced by this change in slope are given by 


P A (Re II) 

./ — 7 r 

p A (ReH) 

J A(ReII) cos 0 dd 


de 

sin 0 dd 
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where A (ReH) is determined from equation (26). These 
changes in slope are effective in eliminating or reducing 
f Re H(e M ) dd and fBe H(e iB ) cos Q dd. (See figs. 6 and 7.) 

Another method for adjusting Re H(e i0 ) is to use a multi- 
plier on the entire velocity distribution; that is, 

K(e)=e- a K{6) 

then 

A (ReH)= 2 cogh 2Jc _ 2 cos 26 

and 

£ 

A (Re If) sin 6 d8= 0 
A (ReH) cos 6 tf0=O 



Changes in the sine and cosine integrals may be obtained 
by using different multipliers over the upper and lower 
surfaces : 


K(8) = K(6)e-°i 
K(6)—K(6) 


9,<9<e n 
d>t<9< 6 1 2tt 


then 


A (ReH) 2 cosh 2k — 2 


cos 26 


0,<0<0„ 


and 


A (Re H) 2 cosh 2 


A ( ReH)dd = 




k— 2 cos 2 6 


e„<e<e l + 2 i t 


di tan 6 

2 sinh 2 k an tanli 


? Y" 

k)o , ' 


( 


& 2 


+ 

j tan 6 Y l+2 * 


2 sinh 2 k tan tanh 1c 


t 


A (Re IT) 


sin 6 dd= — 


__(i\ , , cos_£_V 

4 cosh / ( : an 1 cosh k 


)*■- 


a, , _j cos 6 Vi + 2» 

4 cosh it 1 cosh /rA„ 


.f_/(^^) cos 0 tan_ ' iHTfrX,' 

/a j _i sin 0 V< 

\4 sinh A: an sinh &/ 0 „ 


+ 

l~h27T 


These integrals may be easily evaluated because the inverse 
trigonometric functions all enter into the computation of 
<pi(9) and have been calculated in the determination of k. 

A possibility for making small adjustments in ReH(e ie ) is 
to use a multiplier that is a simple function of 6, such as 


K(6)=K{d) cos 0+c s,n e) 


then 


and 


A (ReH) 


a+b cos 0+c sin 6 
2 cosh 2 k — 2 cos 26 


r A (ReH)d6= 

t/ — 7T 

I A (ReH) sin 6 dd= 

< — 7 T 

p A (ReH) 


cos 6 dd= 


air 

sinh 2k 

ctt (1 — tanh k) 
2 

b-rr (1 — tanh k ) 

2 tanh k 


Various combinations of these methods for adjustment of 
Re H(e ie ) have been used in the illustrative examples of this 
report. 

After the adjustments in Re H(e ie ) have been made, the 
final values should be checked in equations (23), (24), and 
(25), as indicated in table II, lines 42, 43, and 44. At this 
time, K(6) should be computed in order to determine whether 
/v(0)<2, as required by part (d) of equation (5), reference 1. 
If K(9) does not satisfy this inequality, a different, modifi- 
cation of Re II(e iB ) is necessary. 




Figure 7.— Products MR? II) cos 0 and MR? II) sin 0 for k— 0.20 and for A r (0VA'(0)= 0.950. 
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The final values of Re H(e i0 ) and 7j(B) should lie plotted in 
order to insure that no excessive fluctuations have been 
introduced. At this time, q{6) should he compared with the 
original q{6) to make certain that the velocity has not been 
changed hcybnd acceptable limits. Because the potential 
range is fixed, these changes will change, the arc length. If 
desired, 7/(6) may be obtained as a function of the new arc 
length s by 

5(0)=3[0(s)]==5(s) 

and 

s(0) = dd d l <0<8' l + 2ir 

Jo, q \0) 


The values of “-vuw 
<Z W 


arc given in table II, line 50. 


Computation of ImH{e ie ). — When ReH(e ie ) satisfying all 
requirements has been obtained, the computation of the 
conjugate function ImII(e ie ) is done by direct numerical 
integration of Poisson’s integral 


/w/7(e lS )=^- r ReU(c :r ) cot T ~° dr 
J-T 2. 


on the 80 basic points using an extension of the method 
developed by Dr. Glenn H. Peebles. The integration is 
accomplished by replacing ReII(e t0 ) by 40 parabolas on the 
80 basic, points and integrating analytically the product of 
the parabolas and the cotangent term. The two sets of 
80 coefficients so obtained for the integration are given in 
tables VII and VIII. (See appendix A for derivation.) 
The set in table VII is used for the computation of tiie con- 
jugate at the end points of the parabolas; that is, at the even- 
numbered points 0, 2, 4, . . . 78. The set in table VIII 
is used for the computation of the conjugate at the odd- 
numbered points 1 , 3, 5, ... 79. Thus, to obtain ImJI(e ie ) 
at one of the basic points, for example, point 7, the value of 
Refl(e if ) at this point is multiplied by the first coefficient 
(0.000000) in table VIII, the next value of ReH{e ie ) (at 
point 8) is multiplied by the next coefficient (0.412308), 
and so forth, and the sum of these 80 products is the desired 
value of Im H(e ie ) at that point. Hence to obtain I in H(e i0 ) 
at the 80 points requires 80 such accumulative multiplica- 
tions. This computation is done very efficiently on an elec- 
tronic calculating punch using only 415 cards and taking 
approximately 2 hours, including the time for key punching 
and verifying. 

When Re II(e i/r ) has large fluctuations so that some of the 
80 basic intervals arc divided by half-points, the preceding 
values of ImH{e i0 ) must be corrected to take into account 
the difference between the values of ReH(e ts ) and the basic 
parabolas at these half-points. The correction to be added 
to the calculation is the harmonic conjugate of these differ- 
ences. The coefficients for computing the value of the con- 
jugate (at the 80 basic points) of these differences are given 
in table IX. (Sec appendix B.) This computation has been 
arranged on two concentric disks — the coefficients are carried 
by the outer disk and the differences (denoted by 5) are 
entered on the inner disk in the proper places (fig. 8). It is 



easy to see where these differences will make, a significant 
contribution to hnll{e ,e ) and to obtain the contribution by 
accumulative multiplication. This correction is then added 
to the calculated results. 

COMPUTATION OF BLADE COORDINATES 

After the corrected values of the conjugate have been 
obtained, the blade coordinates arc given by equation (17), 
which for convenience of computation may be written in the 
parametric form 



cos [ 0 + ^ +g 2 (0) — tan 1 (27) 

sin [e +^+52 (0) — tan -1 dd ( 2S) 


where 


gi(0) = n ^+/wiC (<?'*) + ( 2 cosh 2 A — 

2 cos 28)ImII{e ie ) 


Those integrals are evaluated by Simpson’s rule. In order 
to obtain sufficient accuracy, use of at least the same number 
of points as were used in the evaluation of Re II(e ie ) is 
advised. Because the values for Im II{e i0 ) were calculated 
at the 80 basic points only, interpolation of g 2 (8) is necessary 
when half-points arc used. The values of x(8) and y(8) are 
given in table II, lines 70 and 71. 

EXAMPLES 

Several examples have been computed to illustrate some 
of the variety of conditions to which the method may be 
readily applied. No attempt has been made to pick the 
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best velocity distribution or to obtain the most desirable blade 
shape, as to do so would lead to many considerations far 
beyond the scope of the present investigation. 

In these examples, the adjustment of the prescribed veloc- 
ity altered the arc length slightly. For uniformity, figures 
9 to 12, which show the final velocity distribution and blade 
shape, were therefore scaled to give an arc length of 2?r. 

Example 1.- — For this example of a cascade blade with 
low turning and a cusped trailing edge («=0), the free-stream 
conditions were taken as 

tq = 0.488 m = 170° 

m 2 =0.478 a*— 180° 

The velocity distribution prescribed on the blade was the 
velocity distribution for an isolated Joukowski airfoil. Ad- 
justments to Re H(e ie ) altered this distribution somewhat 
and the resulting blade (which is very similar to a Joukowski 
airfoil) and velocity distribution are shown in figure 9. 

Example 2. — For an example of an impulse-type blade 
with rounded leading and trailing edges (c=x), the free- 
stream velocities were assigned the values 

■u,=0. 555 ai = 135° 

Uo=0. 555 a, = 225° 

and the prescribed velocity distribution on the blade was 
constant over most of the upper and lower surfaces with 
different values on the two surfaces and varied linearly 




Figure 9. — Final velocity distribution and blade shape /or Example 1. M\ =0.500; «=0 
(cusped trailing edge). 


through the stagnation points. In this example, the expres- 
sion for r (equation (4)) is indeterminant and the value of r 
was obtained by taking the limit as a 2 approaches 225°. 
In this example, complete adjustment of Re H(e ie ) was made 
by multiplying the velocities on the upper and lower surfaces 
by constants. The resulting blade shape and velocity dis- 
tribution are shown in figure 10. A smaller nose or tail 
radius of curvature could be obtained by increasing the slope 
of the velocity distribution through the nose or the tail, 
respectively. 

Example 3. — The free-stream velocities chosen for this 
example were representative of a compressor stage with 

=0.583 a, = 135° 

u 2 = 0.417 a 2 =155° 

In order to keep the changes in the prescribed velocity to a 
minimum, three ratios of lower surface length to upper surface 
length of 0.95, 0.90, and 0.85 were used with the same pre- 
scribed velocity distribution (but, of course, different circle 
flows) and the integrals in equations (23), (24), and (25) were 
quickly approximated using only 20 points. By use of these 
results, a ratio that ■would give the smallest integrals was 
chosen (a ratio of 0.93) and the blade shape with rounded 
trailing edge (e=x) was computed in the usual manner. 



/.Or 


u .5T 


L I i 1 

0 n 2n 

s 

Figure 10.— Final velocity distribution and blado shape lor Example 2. Mi =0.573; e=jr 
(rounded trailing edge). 
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The blade shape and the velocity distribution are shown in 
figure 11 (a). By use of the same circle flow and a change 


in 0,(6), a blade having a 10° 


trailing-edge angle ^ 6== Y§) and 


similar velocity distribution was obtained (fig. 11 (b)). 

Example 4. — In this example of a highly loaded cascade 
blade' with a rounded trailing edge («=7r), the free-stream 
velocities were 


^ = 0.579 «!= 135° 

u 2 = 0.369 a 2 =180 o 


The prescribed velocity distribution, which was similar to 
the distribution on a Griffith airfoil, had an abrupt decrease 
in value on the upper surface for use with suction. The 
blade and velocity distribution are shown in figure 12 (a). 
By use of the same circle flow and a change in 2(0), a blade 
having a cusped tail (e=0) and essentially the same velocity 
distribution was obtained (fig. 12 (b)). 

DISCUSSION 


Re H(e t0 ) will be infinite at 0=0,. Fairing Re H(e iB ) smoothly 
through 6=6, and keeping the values finite will, however, 
insure obtaining the specified angle at the trailing edge of the 
blade. The shape of the blade in the immediate vicinity of 
this point will depend on the velocity prescribed in this 
vicinity. If desired, the blade shape at the tail can easily 
be modified by changing Im H(e ie ), which changes the angle 
of the tangent to the blade (the angle in equations (27) and 
(28)). The corresponding change in Re H(e is ) is computed 
from Poisson’s integral with the constant term zero, and 
the modification in the velocity is obtained from this change 
in Re II(e ie ). The change in Im H(e ie ), denoted by A (Im H), 
should be chosen to satisfy the following conditions: 


j* A(ImIT)dd=0 
|* A(Im H)cosd dd=0 
j A(ImH)sm0 cl9=0 


Specification of the trailing-edge angle of the blade requires 
that u(,s) have a zero of order ~ at s=0 and s = 2ir. If the 

Z7T € 

prescribed velocity does not go to zero in the, proper manner, 


V 



For a prescribed velocity distribution requiring only mod- 
erate changes, the blade shape can be obtained in approx- 
imately 50 computing hours using 80 points and 4-decimal 
accuracy. The. time depends, however, on the degree of 
familiarity with the method and on the extent of the permis- 




Figure 11. — Final velocity distribution and blade shape for Example 3. Jt/i =0.604. 
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(a) «=t r (rounded trailing edge). 

Figure 12.— Final velocity distribution and 

siblc modification to the velocity. When several examples 
have the same free-stream conditions and cascade spacing, 
the time is considerably reduced because the circle flow need 
be computed only once. 

With the calculation based on a spacing of SO points with 
half-points around the nose and the tail, the method gives 
accurate results in all cases in which the parameter k is not 
less than 0.10. Experience has shown that for cascades of 
moderate stagger and turning k will usually be greater than 
0.10 when the solidity is less than 1. In examples 1, 2, 3, 


1.0 r 



O n 2n 


s 

(b) e = 0 (cuspcd trailing edge), 
le shape for Example 4. .\/i = 0.599. 

and 4, the values of k were 0.2G00, 0.1 100, 0.2008, and 0.2851 
and the solidities were 0.91, 1.01, 0.92, and 0.74, respectively. 
In applying the method to other cases, a finer point spacing, 
which would require new coefficients for the integrations, 
should be used. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, August 28, 1950 



DETAILED COMPUTATIONAL PROCEDURE FOR DESIGN OF CASCADE BLADES 

APPENDIX A 

COEFFICIENTS FOR DETERMINING CONJUGATE FUNCTION 


The coefficients for determining Im II (e iB ) are obtained 
by considering the contributions of one of the parabolas to 
the conjugate at one of the basic points. The contribution 
I of the (A + l)th parabola through the points Re H(e ie at), 
Re H{e ie 2 k+ 1 ), and Re II (e ,e 2 k+ 2 ) to the conjugate at any 
point 8= d m is given by 


_ 1 r e 2k+ 2 

-^ == 27r jo ®2t+l) 2 +6(6 — #2*+l) + 

1 ' 

f) 

— 

ReH {e i0 ™+ 0] cot ^ g "‘ c/5 

(2i+9) 

9 + 

- 


where 


(Al) 


lie II(e + 2 ) — 2 Re H(e i0 ™ + >) + Re //(+%) 
a- _ 

, Re II(e '»** + 2 ) - Re //(+%) 

2 A" 

and A is the length of the basic interval. By substituting 
8-2k+\ — 8 ,h -\-lh 

in equation (Al) and expanding in scries, equation (Al) can 
be integrated to give 

I — Ri Re Hie 1 )+S,ReII(e' < "'+«». , ) + 

T, Re Hie 1 *+»-» ) 

where (for 1^ 0, 1) 

2 (F,Z a + VJJ+V.L'^- • • .) + 4 (IF,Z+ 

w'*i*+ir,L*+ ■■■)] 


s,=- 




£ + 1 , 
sin —— A 

In ^£y 8 (W,L+ W 2 L 3 +WJJ-\ ) 

sin ~~ h 


T =if _r 1 MY . 2 MY , 17 MY , 1 

1 A [_ 3 \ 2 / ~^ 15 \ 2 / ‘ 105 \ 2 / 4 ~ 2835 \ 2 / 4 J 
2 (VJJ+V 2 L^+V i L»+ ■ ■ -) + 4 (JF.Z + 

TF 2 Z 3 +PF 3 Z 3 + • • 0 } 


+ 


and 


. A 
tan ^ 

V Hi 

tan ~zr~ 


1+— ' 


tan 2 

v A / 


1 ( tail £) 


L (2.7 + 1) ' (2j+3) 3 

1 ( tan f) 4 , 1 ( tan !)° 


( 2.7 + 5 ) 


(2.7 + 7) 


tan 


W> 


2(tan|)‘ 


L(2.7+ 1) (2j + 3) 


1 

2 .3 (tan §)“ 

1 

44(tan 

(2.7 + 5) 

45 

(2./+7) 105 

1 

56:i(tan|y 



(2.7 + 9) 

1 57 5 




AY 


+ 


S 0 — 0 


When 1=0, 

r.= [1 -i (0'_ (|y_ 


42, 


When /= 1 , 

«-i(r 


A- A 1 


A 0 


'36 1350 26,4(30 425,250 


-•) 


J \=—~ (y 
7T \J 


A 2 2A 4 


8" r 675 _r 8820^ _ 2 


12,625' ) 


11 


'] 


and R 1 is undefined because the integrand becomes infinite. 

The final integration coefficients, denoted by N 0 , A+ 
A+ . . . are obtained from the preceding values of R h S h 
and T h When 9 m is an end point of a parabola (m is even), 
then 


AG=0 

N 2 =Ti j tR 3 


y (A 2 ) 


N 2j—i — S 2 j-i 

A+? = T 2 j_ 1 + R ij+ ! 


where N 0 is determined by Cauchy’s principal value. When 
6 m is a midpoint of a parabola (m is odd), then 
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N 0 =S 0 = 0 
N^To+R* 

n 2 =s 2 

N 2 j-i=T 2 j- 2 -\-R 2 j 

N 2 j=S 2 j 


^ (A3) 


Tlie coefficients for 80 points given in tables VII and VIII 
were obtained from equations (A2) and (A3) using A= 27+80. 
The coefficients for other point spacings can be obtained from 
equations (A2) and (A3) by using the proper value of h. 
Because of the symmetry of the coefficients, only half of 
them need be calculated. 


APPENDIX B 

COEFFICIENTS FOR CORRECTION OF CONJUGATE 


The coefficients used in the wheel correction are obtained 
by determining the contribution of one of the half-point 
differences to the conjugate at one of the basic points. 
Denoting the contribution of the difference &+ ■ (the dif- 
ference at 0 — dk+i) to the conjugate at 0=9 m by /, then 



which, on substituting 


Consequently, the coefficients for the wheel correction 
denoted by Nb, Ab+|, Ab+i, . . ., are 


Ab 


8 (\ h 2 


-§( 


h 6 


AT, + ,=; 


4 288 43,200 3,386,880 

h 


sin (/+!)■ 


In 


. Lh 
sm 2 


-2(E 1 G+E 2 G 3 +E 3 G°+ . . .) 


0*H 


i + 


becomes 


1= 


2 tt Je m +Ih 


4 8, 


n + l 4- 2 


s — (\« + 


The coefficients in table IX, which arc for 80 points, were 
obtained from equation (B2) using h = 27+80. By use of 
other values of h in equation (B2), the coefficients for dif- 
ferent point spacings can be obtained. Because of the sym- 
metry of the coefficients, only half of them need be computed. 

REFERENCE 




+ ,}cot e ~^de 


(Bl) 


1. Costello, George R. : Method of Designing Cascade Blades with 
Prescribed Velocity Distribution in Compressible Potential Flows. 
XACA Rep. 978, 1950. (Formerly NACA TN 1970.) 


Integration of equation (Bl) by series expansions gives 
(for ( + 0) 


1 = 5 


m + 1 + 2 


sin (1+1) ~ 

In — - ~-2(E,G+E 2 G 3 +E 3 G*+ . 


sm 


lh 


- •) 


where 


Q-- 


i h 

tan - 


tan (21+1) - 


and 


( 4 tan 

1 21 

, i 231 


V h ) 

L(2j+1) 3(2.7 + 3) 1 (2j + 5) 

45 


1 441 

( tan i )' 

L 1 

563 ^tan 


(2.7 + 7) 

105 

+2j + 9) 

1 575 

- 

Wlten 1 = 0, 

'-<■*» I < 

n h 2 

/t 4 

h 6 

• .^1 

A 288 

43,200 

3,386,880 

) 


TABLE I— OUTLINE FOR COMPUTATION OF k 


Lino 

i 

Operation 

Remarks 

1 

Assume a value of k 


2 

sinh (1) 


i 3 

cosh (1) 


i 

tanh (1) 


5 

Im A-Ivi B 
2 Re A X W 


6 

tan -1 (5) 

(0) is X; express in degrees; — 90°<X<90° 

7 

see (6) 


S 

i 9 

(Im A Aim B) (2) 


2 Re A ' N (7) 

sin -1 (8) 

(9) is X+0„; -9O°<\+0„<9O° 

10 

(9) -(6) 

(io) is e n 

i ii 

-(10) -2(6) -180° 

(11) is d t 

1 12 

(11) +360° 


13 

sin (10) 


14 

cos ( 10 ) 


! 15 

tan (10) 


! 1(3 

sin (12) 


1 17 

COS (12) 


18 

tan (12) 


1 19 

(18) -5- (2) 

, 

20 

(13)-*- (2) 


21 

(17) +(3) 


22 

(14)-!- (3) 


23 

(18)-H4) 

1 

24 

(lS) + (4) 


25 

tan -1 (19) 

Express in radians; — ^<(25)<-^ 

20 

tan-' (20) 

Express in radians; — ^<(20)<-^ 

27 

tanh -1 (21) 


28 

tanh -1 ( 22 ) 


29 

tan -1 (23) 

Express in radians; same quadrant and 

30 

tail- 1 (24) 

sign as QtA 2tt 

Express in radians; same quadrant and 

31 

(ImA-Im K)X[(25)-(20)J 

sign as 0„ 

32 

—2Re .'IX! (27) — (28)] 


33 

(Im A Aim B)X[(29)-(30)1 

i 

34 

(31)+(32)+(33) 

Must equal <p c (2Tr)—<pc(s„) for correct . 



value of k 




I I 


imm ■■ ii 


I II M 
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Line 


2 

3 

4 

5 

6 
7 

“ 8 " 

9 

10 
11 
12 

13 

14 

15 


17 

18 


19 


20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 


39 


Operation 


Assign values of 0 at 4.5° intervals 
beginning at 0=— *• 

2 ( 1 ) 

COS (2) 
cosh 2fc— (3) 

(4i?e A cosh k sec X) -5- (4) 

(I) +X 

sin, (6). - 

(7) — sin. (0n+X) 

(8) X(5) 
sin (1) 
cos (l) 
tan (1) 

(10) -i-sinh k 

(II) -5-cosh k 
(12)-5-tanh k 

tan-* (13) 

tanh~* (14) 
tan~> (15) 

(Im A—Im B)X(16) — 

2Re AX(17)+ 

(Im A+Iin B)X(18)4-Z* 


Q(fi) 

l+(20) 2 

1+V(21) 

2 (20) + (22) 

9.-0) 
cos (24) 

2-2 (25) 

(26) «/* 

(27)X(23) 

2 w 

V( 29) 

I (9) |X (30) + (28) 

In (31) 

(Cs— Cfl)e i X(10) 
(CV-C0e*X(n) 
[(Gi+Ci)+(33)+(34)]+2 
(32)]-(35)+i]+(29) 


(30) ile 


r 

— n 

J* (36) sii 

— 7T 

J* (36) cos 0 do. 


) sin 6 do 


Remarks 


Half-points may be added later to these 
basic 80 points where necessary 


sec X is line 7, table I 
X is line 6, table I 

sin (0„- f X) is lino 8, table 
(9) is v(0) 


Express in radians;— y< (10) <y 

Express in radians; same quadrant and 
sign as 0 
(19) is <t>i(0) and 

D-~(Im A-Im B) tan-i 4^—;+ 


2 Re A tanh~* 


sinh k 

c os 0t 
cosh k 

(Im A-\-hn B) tan-* tar ) 

* tanh k 

Obtain by matching potentials 

(23) is the distorted velocity IC(d) 
Omit these lines when n = 0 

(28) = (23) when «=0 


(35) is Re C(e i8 ) 

(36) is Re H(e'*) 

Evaluate by Simpson's rule 


Evaluate by coefficients in tables III and 
V 


Evaluate by coefficients in tables IV and 
VI 


Line 


41 

42 


44 


45 

46 

47 

48 

49 


51 

52 

53 


55 


58 

59 

60 
61 
62 

63 

64 


65 

66 

67 

68 
69 


Operation 


A(Re H) 
(36)+(40) 
f r (ii) <te 


(41) sin 0 do 


(41) cos 0 do 


— (40) X (29) 

(46)X(23) 

4- (47)2 
4 (47) 4- (48) 

I (9)1 +(49) 

(Cs-C,)e‘X( 11) 
(Cj-C,)£*X(10) 
l(C,+C«)+(51)- j (52))+2 
Uncorrected Im H(e is ) 

Correction for Im H(e ia ) 


(54)+(.W) 

|(24)+180°]x| 

(29)X(56) 

57.29578 [(58) +(53)] 
(57) + (59) 
sin (2) 

(3)-«* k 

(61) - 5 - (62) 
tan-* (63) 


(1) +(60) +90° -(64) 
cos (65) 
sin (65) 

(50) X (66) 

(50) X (67) 


r 

f°m 


(68) do 


do 


Remarks 


Changes in Re H(e* B ) to make integrals 
zero; see text 

(41) is Re H(e“) 

Must be zero 

Must be zero 

Must be zero 


(47) is K(0 ) ; must have K(0) <2 
(equation (5d), reference 1) 

(49) is q (0) 

(5°) is | 


(53) is lm C(e ■») 

Obtained from machine calculation using 
coefficients in tables VIII and IX 

Obtained from differences between Re //(«•*) 
and parabolas using coefficients from 
table IX 

(56) is Im IJ(e iB ) 

(57) =0 when «=0 


(60) is qi (0) in degrees 


Express as a positive angle in degrees; 
quadrant chosen so that sin (64) and cos 
(64) have same signs as (61) and (62), 
respectively 


(70) isa:( 0 ); evaluate integral by Simpson’s 
rule 


(71) is y(0 ); evaluate integral by Simpson’s 
rule 


TABLE III — WHOLE-POINT COEFFICIENTS FOR SINE 
INTEGRAL 



0 

10 

20 

30 

40 

50 

60 

70 

0 

0.00000 

-0. 03707 

-0. 05242 

-0. 03707 

0.00000 

0. 03707 

0. 05242 

0. 03707 

1 

— .00821 

-.07958 

— . 10433 

-.06797 

. 00821 

. 07958 

. 10433 

. 06797 

2 

-.00820 

-.04241 

-.05178 

-.03081 

.00820 

.04241 

. 05178 

.03081 

3 

-.02443 

-.08923 

-.10176 

-. 05468 

.02443 

. 08923 

. 1017G 

. 05468 

4 

-.01620 

-. 04671 

-.04986 

-.02380 

.01620 

. 04671 

. 04986 

.02380 

5 

-.04005 

-.09669 

-.09669 

-.04005 

.04005 

.09669 

.09669 

.04005 

6 

-.02380 

-.04986 

-.04671 

-.01620 

.02380 

. 04986 

. 04671 

.01620 

7 

05468 

-. 10176 

-. 08923 

-.02443 

. 05468 

. 1017C 

.08923 

. 02443 

8 

03081 

-.05178 

-.04241 

—.00820 

. 03081 

. 05178 

.04241 

.00820 

9 

-.06797 

10433 

— . 07958 

-.00821 

. 06797 

. 10433 

. 07958 

.00821 


TABLE V— HALF-POINT COEFFICIENTS FOR SINE 
INTEGRAL 



0 

10 

20 

30 

40 

50 

60 

70 

0 

0. 00000 

-0. 01852 

-0. 02619 

-0. 01852 

0.00000 

0. 01852 

0. 02619 

0. 01852 

H 

-. 00206 

-. 03844 

~. 05231 

-.03554 

.00206 

. 03844 

.05231 

. 03554 

1 

— . 00206 

-. 01991 

-. 02611 

01701 

.00206 

.01991 

.02611 

. 01701 

1J4 

— . 00615 

-.04111 

— . 05199 

— . 03241 

.00615 

.04111 

.05199 

. 03241 

2 

— . 00410 

— . 02119 

-. 02587 

—.01539 

.00410 

. 02119 

.02587 

.01539 

2 a 

01021 

-. 043.53 

-. 05135 

-.02908 

. 01021 

. 04353 

. 05135 

.02908 

3 

-. 00611 

02233 

02547 

01368 

. 00611 

.02233 

.02547 

. 01368 

3 Vz 

-.01421 

— . 04568 

— . 05039 

— . 02558 

. 01421 

.04568 

. 05039 

.02558 

4 

00809 

02333 

-.02491 

-.01189 

.00809 

.02333 

. 02491 

. 01189 

i'A 

-.01812 

04754 

-. 04912 

-. 02192 

. 01812 

. 04754 

. 04912 

. 02192 

5 

-. 01002 

-. 02419 

-. 02419 

-. 01002 

. 01002 

. 02419 

. 02419 

. 01002 

5'A 

-. 02192 

— . 04912 

-. 04754 

01812 

. 02192 

. 04912 

. 0475'! 

.01812 

6 

-. 01189 

-. 02491 

02333 

-. 00809 

.01189 

. 02491 

. 02333 

. 00809 

BA 

-. 02558 

-. 05039 

-. 04568 

01421 

.02558 

.05039 

.04568 

. 01421 

7 

-.01368 

-. 02547 

-. 02233 

-. 00611 

.01368 

.02547 

.02233 

.00611 

7 H 

— . 02908 

-. 05135 

-. 04353 

-. 01021 

.02908 

. 05135 

.04353 

. 01021 

8 

— . 01539 

-.02587 

-.02119 

-. 00410 

.01539 

.02587 

.02119 

.00410 

8'A 

-.03241 

-. 05199 

-.04111 

-.00615 

.03241 

.05199 

.04111 

.00615 

9 

-. 01701 

-.02611 

01991 

-.00206 

.01701 

.02611 

.01991 

.00206 

9 H 

— . 03554 

05231 

03844 

— . 00206 

.03554 

.05231 

.03844 

.00206 


TABLE IV— WHOLE-POINT COEFFICIENTS FOR COSINE 

INTEGRAL 



0 

10 

20 

30 

40 

50 

60 

70 

0 

-0. 05242 

-0. 03707 

0. 00000 

0. 03707 

0. 05242 

0. 03707 

0. 00000 

-0. 03707 

1 

-.10433 

00797 

.00821 

.07958 

. 10433 

. 06797 

-.00821 

-.07958 

2 

-.05178 

-.03081 

.00820 

. 04241 

. 05178 

.03081 

-.00820 

-.04241 

3 

-. 10176 

-.05468 

. 02443 

. 08923 

. 1017(5 

. 05468 

02443 

-.08923 

4 

-.04986 

-.02380 

. 01620 

. 04671 

. 04986 

.02380 

-.01620 

-.04671 

5 

09669 

-.04005 

. 04005 

. 096(59 

. 096(59 

. 04005 

-.04005 

-.09669 

6 

-.04671 

-.01620 

.02380 

. 04986 

. 04671 

. 01620 

-. 02380 

-.04986 

7 

-.08923 

-.02443 

. 054(58 

. 10176 

.08923 

. 02443 

-.05468 

-.10176 

8 

-.04241 

-.00820 

. 03081 

.05178 

. 04241 

.00820 

-. 03081 

-.05178 

9 

-. 07958 

-.00821 

. 06797 

. 10433 

. 07958 

.00821 

-.0(5797 

-. 10433 


TABLE VI— HALF-POINT COEFFICIENTS FOR COSINE 

INTEGRAL 



0 

10 

20 

30 

40 

50 

60 

70 

0 

-0. 02619 

-0.01852 

0. 00000 

0. 01852 

0. 02619 

0. 01852 

0. 00000 

-0. 01852 

'A 

05231 

-. 03554 

.00206 

.03844 

.05231 

. 03554 

00206 

03844 

1 

-.02611 

-. 01701 

. 00206 

.01991 

.02611 

. 01701 

-. 00206 

-.01991 

i M 

—.05199 

-. 03241 

.00615 

.04111 

. Q5199 

. 03241 

— . 00615 

-.04111 

2 

-. 02587 

-. 01539 

. 00410 

. 02119 

.02587 

. 01539 

00410 

-.02119 

2 'A 

05135 

-.02908 

. 01021 

. 04353 

. 05135 

.02908 

-.01021 

-. 04353 

3 

-.02547 

-. 01368 

.00611 

.02233 

.02547 

. 01368 

00611 

-. 02233 

m 

-. 05039 

-. 02558 

. 01421 

.04568 

. 05039 

.02558 

-. 01421 

-. 04568 

4 

-. 02491 

-.01189 

. 00809 

.02333 

. 02491 

.01189 

00809 

-. 02333 

4>2 

-.04912 

-. 02192 

.01812 

.04754 

.04912 

. 02192 

-.01812 

-. 04754 

5 

-. 02419 

-. 01002 

. 01002 

. 02419 

. 02419 

. 01002 

-. 01002 

-. 02419 

5'A 

04754 

-. 01812 

. 02192 

. 04912 

. 04754 

.01812 

-.02192 

-.04912 

6 

-.02333 

00809 

. 01189 

. 02491 

.02333 

. 00809 

-.01189 

— . 02491 

6 H 

-. 04568 

01421 

.02558 

. 05039 

.04568 

. 01421 

02558 

05039 

7 

-.02233 

-. 00611 

.01368 

.02547 

.02233 

. 00611 

-. 01368 

02547 

I'A 

-.04353 

— . 01021 

.02908 

. 05135 

.04353 

- 01021 

— . 02908 

-.05135 

8 

02119 

-.00410 

.01539 

.02587 

.02119 

. 00410 

-.01539 

02587 

8'A 

-.04111 

-. 00615 

.03241 

.05199 

.04111 

.00615 

-.03241 

-.05199 

9 

-. 01991 

-. 00206 

. 01701 

.02611 

.01991 

.00206 

-. 01701 

02611 

9'A 

-.03844 

. -00206 

.03554 

.05231 

.03844 

.00206 

-. 03554 

— . 05231 
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TABLE VII— COEFFICIENTS FOR CONJUGATE AT EVEN- 
NUMBERED POINTS 



0 

10 

20 

30 

40 

50 

60 

70 

0 

0. 000000 

0. 020031 

0. 008323 

0. 003449 

0. 000000 

—0. 003449 

—0. 008323 

—0. 020031 

1 

. 636402 

. 036217 

.015415 

. 006151 

-. 000655 

-.007686 

018042 

-.045294 

2 

. 050356 

. 016305 

.007110 

. 002706 

000655 

004243 

-. 009743 

-. 025472 

3 

. 144119 

. 029799 

. 013146 

. 004702 

-.001973 

-.009338 

-.021159 

059345 

4 

. 050944 

.013568 

. 006049 

. 001999 

— . 001319 

-.005102 

-.011449 

-. 034277 

5 

. 084480 

. 024968 

.011141 

. 003316 

-. 003316 

011141 

-. 024968 

-. 084480 

6 

. 034277 

.011449 

.005102 

.001319 

-.001999 

-.006049 

-.013568 

-.050944 

7 

. 059345 

.021159 

. 009338 

.001973 

-.004702 

-. 013146 

-. 029799 

-. 144119 

8 

. 025472 

. 009743 

. 004243 

. 000655 

-. 002706 

-.007110 

-.016305 

— . 050356 

9 

. 045294 

. 018042 

. 007686 

. 000655 

-.006151 

-.015415 

-.036217 

-.636402 


TABLE VIII— COEFFICIENTS FOR CONJUGATE AT ODD- 
NUMBERED POINTS 



0 

10 

20 

30 

40 

50 

60 

70 

0 

0. 000000 

0. 040322 

0. 016677 

0. 000906 

0. 000000 

-0. 006900 

—0. 016677 

—0. 040322 

1 

. 412368 

.018011 

. 007694 

. 003072 

000327 

-. 003839 

-. 009002 

022467 

2 

, 223706 

. 032759 

. 014242 

.005417 

-.001312 

008495 

019528 

-. 051462 

3 

. 065520 

. 014841 

.006563 

. 002349 

-. 000986 

-.004663 

-.010554 

— . 029282 

4 

. 106592 

. 027228 

.012115 

. 004003 

— . 002641 

010218 

-.022960 

-. 089820 

5 

.041111 

. 012446 

. 005563 

. 001657 

001657 

-. 005563 

-. 012446 

-.041111 

6 

. 059820 

. 022960 

. 010218 

. 002641 

-.004003 

-.012115 

-. 027228 

-. 106592 


. 029282 

. 010554 

. 004663 

. 000986 

-. 002349 

-. 006563 

-. 014841 

-.065520 

S 

. 051462 

.019528 

. 008495 

. 001312 

-. 00.541 7 

-. 03 4242 

032759 

-. 223706 

9 

. 022467 

. 009002 

. 003839 

. 000327 

-. 003072 

-. 007694 

-.018011 

-.412368 


TABLE IX— COEFFICIENTS FOR CORRECTION TO CON- 
JUGATE 



0 

10 

20 

30 

| 40 

50 

60 

1 70 

' 0 l 

0. 6.36565 1 

0. 019061 

0. 008014 

0. 003262 

1 

1 -0.000104 

-0.003645 

-0.008669 

-0.021304 

; i 

.11401 1 

.017188 . 

. 007406 

. 002890 

1 -.000492 , 

-.004042 j 

-.009379 

-.024049 

■ o 

1 .085301 

| .015596 

. 006840 

. 002528 

' —.000821 

— . 004455 

—.010156 

027496 

' 3 

.060498 

.014222 

. 006309 

. 002175 

. -.001153 

004885 

| -.011011 

-.031974 

4 

1 .046782 

! .013021 | 

. 005809 

. 001829 

-.001488 | 

-.005.335 ■ 

-.011960 | 

-.038045 , 

1 5 

, .038045 

.011960 1 

. 005335 

.001488 

1 -.001829 

-.005809 i 

-.013021 

-. 046782 

I 6 

.031974 

i .011011 

. 004885 

.001153 

• -.002175 | 

1 -.006309 1 

! -.014222 

-.060498 , 

1 ~ 

1 .027496 

1 .010156 ; 

. 004455 

. 000821 

. -.002528 : 

-.006840 . 

-.015596 

1 -.085301 

8 

.024049 

.009379 ! 

. 004042 

.000492 

-.002890 1 

-.007406 

-.017188 ■ 

-. 144611 

1 9 

.021304 

. 008669 

. 003645 

. 000164 

-.00321)2 

-.008014 1 

1 

-.019061 1 

| -.636565 , 
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